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Abstract. Thermodynamic properties of the helical spins with magnetoelectric 
coupling is simulated by Monte Carlo algorithm. It is shown that the spin-spin 
correlations are affected by the thermal fluctuations and the electric held. Results 
for the temperature-dependence of ferroelectric polarization and the susceptibility are 
also presented in connection with the multiferroic phase transition. 
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1. Introduction 

The efficient control of magnetism in terms of the electric field in the solid has been 
highly desired for many years due to their potential application and fundamental 
curiosity. By now the most straightforward way is to enhance the magnetoelectric 
(ME) coupling in the multiferroics, i.e., materials with both ferroelectric (FE) and 
magnetic orders pp. Here the ferroelectricity results from the loss of inversion symmetry 
induced by the spatial variation of the spin order. Recent experiments indicate that both 
collinear and noncolinear magnetic order can have a potential to produced a spontaneous 
FE polarization P [2] . Among them, the multiferroics with spiral spins are particularly 
interesting [3]. Every nearest-neighbor spin pair produces a local Pj, which can be well 
described by the spin-current model [I], 



with eij being the unit vector connecting two neighboring spins Si and Sj. Since the 
FE behavior is highly sensitive to the applied electric field, the magnetic action is 
then converted to the electric response. The static ME effect has been successfully 
investigated experimentally [5] and theoretically [6] for the helimagnet oxides, which 
would be expected to bring about entirely new concepts in the design of next generation 
devices. On the other hand, the dynamical exchange-striction induces a biquadratic 
interaction between the spins and transverse phonons [7j resulting in the hybridization 
between the spin wave and the fluctuation of the electric polarization, which leads to a 
new collective ME excited mode, so-called electromagnon, providing another promising 
way to enhance the strength of ME coupling by using the resonance in the dynamical 
regime jH |9] . In this paper, the spin-current model is revisited by using the Monte-Carlo 
simulation with classical Heisenberg spins [ID] . We show how the spin-spin correlations 
are affected by the thermal fluctuations and the electric field. It is confirmed that the 
long-range multiferroic order is stabilized by applying electric field that leads to a FE 
phase with spontaneously broken spin chiral symmetry. The polarization susceptibility 
shows an electric field-dependent phase transition. 

2. Spin current model 

We start with a one-dimensional spin chain along the x-axis with a frustrated spin 
interaction. Taking account of electric energy, an effective Ji — J<i model with the ME 
coupling is introduced as follows, 



where E = (0, 0, E) is the external electric field. The notation {ij) nn indicates nearest- 
neighboring (nn) i and j, and (lm) nnn corresponds the next-nearest-neighboring (nnn) 



Pi ~ x (Si x Sj) 




(1) 



(2) 
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/ and m. The competition between the nn ferromagnetic interaction {—J\ < 0) and the 
nnn antiferromagnetic interaction ( J 2 > 0) leads to magnetic frustration and to a spiral 
spin ordering with the pitch angle cos if) = Ji/4J 2 for J 2 > Ji/4 at T = [II] . Without 
external fields, the helical order disappears for T ^ because of the rotational SU(2) 
spin symmetry, only algebraically decaying short-range order is possible. However, the 
ME interaction term, originating from a spin-orbital coupling, is loss of the inversion 
symmetry along the chain. In the presence of an external electric field, the SU(2) spin 
symmetry is broken down to SO(2)xZ2. The spin chirality, 

Ci = Si X Dj+i 

order corresponds to the spontaneous breaking of the discrete Z 2 symmetry and can 
have a nonzero expected value [H]. The electric dipole moments, directly relating to 
the spin chirality, become ordering and we have finite macroscopic FE polarization. 

3. Results and discussion 

The Monte-Carlo simulation is performed on a helimagnet chain with 250 three- 
dimensional classical Heisenberg spins. Spins are randomly updated with the coherent 
rotation, nucleation and Ising-like inversion for each Monte Carlo step in order to 
guarantee that all possible reversal mechanism may occurs in the system and can be 
simulated efficiently [12]. This is especially important in dipolar and chiral system as 
these interactions favors large angles between neighboring spins. Let us now look at 
the data for the static correlation functions (So ■ S n ), n — 1, 2, . . . , 50 for Ji/J 2 = 2, 
cf. FigJU The algebraically decaying of the helical order is clear reflected in the 
amplitude of the correlation function. As expected, the correlation becomes stronger 
as the temperature decreases. The quantity of more fundamental static spin-structure 
factor can be obtained by the Fourier-transformed spin correlation function: 

M 

(S k ■ S-k) = 1 + y~](So ■ S n ) cos kn (3) 

n=l 

Theoretically, M should be infinite but in actual calculations we have taken M = 50. 
The pitch angle is located where the static spin-structure factor is peaked. Figj2] shows 
the maximum of {Sk ■ S_k) tends to diverge at T = 0, indicating the helical long-range 
order. At low temperature, the ground state k = if) and a few around excited states are 
accessible. As temperature increases, the contribution of the critical wave vector k — if) 
decreases. The correlation length of the chirality order shrinks as well, we thus have 
more multiferroic domains. On the other hand, applying even a small electric field can 
elongate the FE order (Figj3]). Unfortunately, the long-range FE order will destroyed 
by the thermal fluctuations too. The transition temperature is electric field-dependence 
[13J. In FigJU we present the result of the polarization (spin chirality) susceptibility. 
For E = O.IJxj the susceptibility shows a broad maximum at T = 0.25Ji. 
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4. Conclusion 

In summary, the thermodynamics of multiferroic spin chain is studied by Monte-Carlo 
simulation based on the classical J\ — J 2 model. We show that the low temperature 
excitations are dominated by the multiferroic domain wall. The ME polarizing procedure 
can be applied to lift the spin chiral symmetry, and then produce a single multiferroic 
domain. The role of the exchange-striction, magnetic anisotropy, and FE dipolar 
coupling for the ME thermodynamics will be investigated further. 
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Figure 1. Spin-spin correlation function as a function of the distance and the 
temperature in the absence of electric field. The angle between neighboring spins 
is ijj = 7r/3 given by J1/J2 = 2. The solid lines are guides for the eye only. 




Figure 2. The static spin-structure factor (Sk ■ S-k) for different temperatures Tj J\ 
= 0.0, 0.1, 0.2, 0.3, 0.5 without electric field E. The pitch angle is demonstrated by 
the maximum at k = tt/3. 
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Figure 4. The total FE polarization (black) and susceptibility (blue) versus 
temperature for the electric field E — O.lJi. 



